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Abstract 



The six-quark instanton induced 't Hooft interaction, which breaks the 
unwanted Ua(1) symmetry of QCD, is also a sourse of semi-classical correc- 
tions to the low energy effective action. It is argued that there emerges a 
dimensionless expansion parameter that introduces a new mass scale, A| ~ 
6 GeV 2 , in the 0~,0 + channels. This scale plays a similar role as the large 
critical mass M c 2 rit > 4.2^-6.6 GeV 2 discovered in the framework of QCD sum 
rules in the 0~, + gluonic channels. In particular, it allows to resolve the rf 
puzzle. To extract A E we calculate the masses of the lightest pseudoscalar 
meson nonet by using the Nambu - Jona-Lasinio (NJL) type Ul(3) x Ur(3) 
chiral symmetric Lagrangian together with the 't Hooft determinant. The 
mechanism which leads to the large value of Ae is scrutinized. 
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1. Introduction 

Unfortunately, there is at present no quantitative framework within QCD to 
deal with its large distance dynamics. The physics of hadrons is approached 
through phenomenological parametrizations usually based on some simple 
ansatz with solid symmetry grounds. There are two important experimental 
facts to support this line of investigation. First, it is known that the chiral 
symmetry of the massless QCD Lagrangian, which should be a good approx- 
imation for light quarks (u,d,s), is not seen in the hadronic spectrum (the 
SU(3)v degenerate multiplets with opposite parity do not exist). It means 
that the hadronic vacuum is not symmetric under the chiral group. Second, 
it is seen from the mesonic spectrum that the £^a(1) symmetry of the QCD 
Lagrangian is badly broken. The SU (3) singlet pseudoscalar 7/ is too heavy 
to be the ninth Goldstone boson. The £7a(1) anomaly is responsable for the 
7/ — n, K, 77 splitting p. It has been understood later [2] that the 1/N C ex- 
pansion can be a relevant approximation to generate hadronic bound states 
and to find the singlet-octet splitting as a next to the leading l/N c order 
effect. 

A qualitatively correct picture of both spontaneous chiral symmetry break- 
ing and Ua(\) breaking at low energies, which is also compatible with the 
conclusions coming from the large N c expansion, is given by instantons E] ■ 
The semi-classical theory based on the QCD instanton vacuum provides con- 
vincing evidences that 2iV/-quark interactions (JV/ is the number of flavours) 
actually exist in QCD and that in the leading l/N c order they are described 
by the 't Hooft determinant jSj, 

C-2N f = «(det qP L q + det qPiiq) (1) 

where the matrices Pl,r — (1 T 7s)/2 are projectors and the determinant is 
over flavour indices. We assume here that all interactions between quarks 
can be taken in the long wavelength limit where they are effectively local 4 . 
At next to the leading 1/N C order this vertex is modified by the tensor term 
which we have omitted in eq.(JT|). Even in this essentially simplified form 
the determinantal interaction contains all necessary features to describe the 
dynamical symmetry breaking of the hadronic vacuum and explicitly breaks 
the axial Ua(1) symmetry. In the following we will asume that quark fields 

4 To lowest order in 1/N C , hadronic physics can be described as a tree approximation 
to some local Lagrangian, with local hadron fields and local interaction vertices [S]. 
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have colour, N c = 3, and flavour, Nf = 3, indices which range over the set 
i = 1, 2, 3. The coupling constant k is a dimensional ([«] = GeV -5 ) negative 
parameter with the large N c asymptotic k ~ 

On lines suggested by multicolour chromodynamics, however, it can be 
argued [2| that the Ua(1) anomaly is negligible in the large N c limit, the 
deviation of the singlet - octet mixing angle from its ideal value is sup- 
pressed, so that mesons come degenerate in mass nonets. Thus, the leading 
order mesonic Lagrangian must inherit the Ul(3) x Ur(3) chiral symmetry 
of massless QCD. To specify the corresponding part of the effective quark 
Lagrangian we consider the four quark NJL type interections [7] 



C 



G 



NJL 



{qKqY + {qilsKqY , (2) 



where A a , a = 0, 1, ... 8 are the normalized (trA a Ab = 2<5 ab ) Gell-Mann matri- 
ces acting in flavour space. The positive four quark coupling G, [G] = GeV -2 , 
counts as G ~ 1/N C and therefore the Lagrangian (|2j) dominates over £q at 
large N c . 

There is another approach to the E/a(1) problem which is also based on 
the method of effective Lagrangians jH]. The way we follow in the present 
work reflects the quark structure of light pseudoscalar and scalar mesons and 
has a built-in mechanism for dynamical chiral symmetry breaking. A similar 
model has been considered for two flavours in [0| and for three flavours in 
[TIH ITTj and has been widely explored since that time ^2] • 

We have used the multicolour asymptotics to motivate our choice of many 
quark interactions (JTJ and (J2J). However, we are not going to follow explicitly 
the idea of 1/N C expansion. As is well known from QCD sum rules, the 
channels with quantum numbers J p = + ,0~ are strongly coupled to the 
gluonic world ^21- It has been argued there that for these channels the 
pictures emerging at N c — >• oo and N c = 3 seem to be qualitatively different 
from each other and the accuracy of the 1/N C expansion becomes worse. 
From a pure phenomenological point of view it would suffice to mention 
here the large deviations from the Zweig rule in the pseudoscalar channel, 
or the well known rj' puzzle: the mass of this meson, being of order 1/N C , 
is unexpectedly too large. QCD sum rules relate these deviations from the 
1/N C counting to the increasing of the mass scale M cr i t characterizing the 
breaking of asymptotic freedom in the corresponding channel, in such a way 
that at large energies E 3> M crit , the 1/N C counting rehabilitates itself. The 
mass scale relevant for rf physics is M crit > 2.0 4- 2.6 GeV jT31 . 
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Since the critical mass is too large, one cannot rely much on the 1/N C 
expansion. Here we use another approximation, namely the stationary phase 
method for the semi-classical path integral bosonization of the effective quark 
Lagrangian 

£ = qii^dft - rh)q + £ N jl + C 6 . (3) 

In this approximation the interactions £njl and £ 6 are considered as contri- 
butions of the same order in h. 

One can also try to improve the lowest order SPA result by taking into ac- 
count the Gaussian fluctuations of the six-quark 't Hooft determinant around 
the stationary phase trajectory [HJ[T3]. If ^6 is not small, corrections can be 
much larger than one could predict starting from 1/N C , and be important for 
the mesonic (0 + ,0~) mass spectra. The dimensionless parameter ensuring 
the smallness of semi-classical corrections in the model is 



c 32G 3 ~ m u 



c 



where Q is the volume of a small Euclidean spacetime box with a side 27t/Ae- 
For certain, we are dealing here with a small effect. However, one cannot 
neglect such quasi-classical contributions until the mass scale Ae associated 
with it is established. The strong 1/N C suppression of (jlj) makes room for 
a large value of still leaving ( small enough. We argue here that the 
cut-off Ae is in a sense similar to the mass scale M crit advocated in ^3] • Our 
assertion is based on a calculation of the mass spectrum of pseudoscalar (0~) 
mesons from which one can extract Ae and show that it is large Ae > 2 GeV. 

It is worth noting that the model under consideration contains a second 
dimensionless parameter [T5] . 

IkIA 1 , \ 

(5) 



AG 2 N c 

with A = m — rh, where m stands for the constituent and rh for the current 
quark mass. The series expansion in e closely corresponds to the 1/N C ex- 
pansion of the model. However, the fit to the meson mass spectrum shows 
that e ~ 0.7, being in contradiction even with the N c = 3 estimate from (jSJ). 
This relatively large (in comparison with 1) value implies large 1/N C correc- 
tions which convert the 1/N C series into a badly convergent one. Sharing 
ideas of paper [TH], we explain this behaviour of the series by the existence 
of a large critical mass in the channel. In this sense, the model perfectly 
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reflects the known resolution of the rf puzzle by generating quasi-classically 
a large mass scale parameter which makes rj light on its natural scale 



2. Characteristic scale of semi-classical corrections 

Let us briefly recall some results of bosonization of the quark Lagrangian 
(J3J). Important details can be found in On the first stage one should 
linearize the many fermion vertices by introducing auxiliary bosonic fields. 
The pure quark Lagrangian, C, is transformed to a mixed meson-quark one 

£mix(<?, 0, a) = C q + C r + A£ r . (6) 

The first term describes the tree level interactions of constituent quarks 
with pseudoscalar, 4> a (x), and scalar, a a (x), U(3) flavour nonets 

C q = qii^dfi -m-a- 175^)9. (7) 

The second term is the leading order stationary phase result 

C r = % tr {U st Ul) + ~ tr ( WUl + W^U st ) 
12 o 

= Ka a + -h ( ^a a a b + -h^Mb + ■■■ ■ (8) 
Here we used the stationary phase condition 

GU a + W a + —A abc Ulul = (9) 

where the totally symmetric constants A abc are defined through the flavour 
determinant det W = A abc W a WbW c . Our notations are the following. The 
trace is taken over flavour indices, any flavour matrix written without open 
index is undestood as summed with the Gell-Mann A a matrices (a = 0, 1...8), 
for instance, W = W a X a . The mesonic fields are grouped in the covariant 
combinations W a = a a + A a — i(j) a . The field U st represents the exact solution 



of the stationary phase condition (JOJ), which we seek in the form U a = s a — ip a 
expanding s a ,p a in increasing powers of bosonic fields a , a a . 

= h a + h { 2a b + h { 2 c a b a c + h^Mc +■■■ (10) 
Pa =h^U b + h { l4> b a c + ... (11) 
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with coefficients explicitly depending on the quark masses and coupling 
constants G, k. The coefficients are fixed by the series of coupled equations 
following from © and obtained by equating to zero the factors before inde- 
pendent combinations of mesonic fields. Due to recurrency of the considered 
equations all coefficients are determined once the first one, h a , has been ob- 
tained [To] . 

An alternative form to the exact solution of eq.flHJ) is provided by the 
1/N C expansion which gives the stationary phase solution, L^", in form of the 
series 

U? = -3 (W a + J^A^WlW} + 0{l/Nl)) , (12) 
yielding for C r 

C r = ~ tr (WW*) - j^- (det W + det W^) + 0(1/N C ). (13) 

In fact, the large- N c limit forces a series expansion for the coefficients of 
Lagrangian (jHJ) and this is how the dimensionless parameter (JSJ reveals itself. 
To see this, let us assume for a moment that SUf(3) flavour symmetry is 
preserved. In this case the only coefficient among h a which is different from 
zero is ho and it is given by 

If the ratio e = kA/(4G 2 ) is small, one can expand the square root inside ho 
in powers of e. This automatically leads to the same expansion for all other 
coefficients in (JSJ) and, finally, one obtains Lagrangian ()13j) at leading order 
in e. This limit is not affected by the semi-classical corrections AC r , because 
these are at most of order ~ 1/N C . One can argue, however, that numerically 
e ~ 0.7 (see our discussion in Sec. 3) what is relatively large, to support the 
fast convergence of the series. 

The third term in © is the next to the leading order correction in the 
semi-classical expansion of the bosonized Lagrangian [To] . 



DC 



71=1 

where 



A£ r = -fT 1 £ K -^J- tr [F aP {(f>, o)} n (15) 



_/,(l)sst /,(l)„st 
ac e ac ire 




st I ' 



(16) 



a/3 
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with s„ = s s ^ — h a . The factor Q 1 may be written as an ultraviolet divergent 
integral regularized by introducing a cut-off A E 5 

/■A E /2 H 4 K. A 4 

One can speculate about the value of the cut-off. By definition Ae belongs 
to the "quark territory" Ae > 1 GeV. On the other side, only effects that 
go beyond standard perturbation theory are included into the Lagrangian 
Cq. Therefore, the characteristic volume of quantum fluctuations, Q, can 
not deviate much from the size determined by non-perturbative fluctuations 
corresponding to classical solutions of the non-linear Yang-Mills equations, 
i.e. instantons. The relevant mass scale is generated by the gluon vacuum 
condensate. Thus, we have crude estimate 

Q ~ (0\—G a G a \0) ~ (330 MeV)~ 4 , A E ^ 2.1 GeV. (18) 

The crucial question is, however, whether this large value is in agreement 
with the general idea of a semi-classical expansion: corrections must be small 
in comparison with the leading order result. At first sight, it seems that we 
have just the opposite case. The factor Q^ 1 ~ A E and one can naively 
expect that large values of Ae will severely break the convergence of the 
quasi-classical series. 

This is not entirely true. To clarify the point and to learn one important 
feature of Lagrangian ()15|) . let us use again the 1/N C expansion. The first 
terms of (JT5j) in this framework are 



AC r -- 



8(2G) 1 



tr (WW r ) + (det W + det W ] ) + C(l/A^ c 2 



(19) 



One can make here a few interesting observations. First, the correction to the 
result ()13|) starts from a term of order. Due to fine cancellations in (|15j) 
it is two orders less than one would naively expect. Second, every term in 
eq. (fT9|) is suppressed by the factor ( (for its definition see eq.(jlj)) as compared 
with corresponding terms in (|13jl. This dimensionless parameter measures 
the size of the semi-classical contribution. Its value must be small, but not 

5 One should not confuse this parameter with the standard ultraviolet NJL cut-off 
A ~ 1 GeV for the quark loops, which represents the mass scale of spontaneous chiral 
symmetry breaking and determines the value of the quark condensate to leading order. 
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necessarily so much suppressed as it follows from the above estimations. 
Actually, nothing forbids us to suppose that such a suppression is partly 
compensated due to the increase of the mass scale Ag, in such a way that £ 
still remains small, ( ~ 10 _1 . 

3. A model estimation of A E 

Let us see how our expectations look numerically. To check our guess one 
should turn to the calculation of the pseudoscalar mass spectrum and ex- 
tract Ae by confronting model results with experimental data. To make our 
consideration not too overfilled with details of numerical calculations we give 
here only final results. We postpone the details until a future publication 

In the following we shall consider the case of SU(2)j x U(l)y symmetry, 
i.e. we take rh u = rhd 7^ rh s . Accordingly, there are alltogether six param- 
eters, rh u , rh s , G, k, A and Ae, shown in Table 1. The last two collumns 
present the dimensionless parameters of the model, ( and e, related to the 
different types of power expansions. In Table 2 are the results of our calcula- 
tions of the pseudoscalar spectrum, together with the weak decay constants 
fn, fx and mixing angle V in the singlet - octet basis (0o,0s)- Inputs are 
indicated by (*). The Latin letter labels on the left hand side identify the 
sets in the tables. 

TABLE 1 . Main parameters of the model given in the following units: [m] = MeV, 
[G] = GeV -2 , [k] = GeV -5 , [A] = GeV. Sets (a, b) correspond to the leading order 
SPA. Sets (c, d) include semi-classical corrections. Constituent quark masses mi 
and corresponding leading order results m,j are also given. 





m u 


m u 




rh s 


O 

m s 




G 


— K 


A 


Ae 


c 


e 


a 


5.3 


315 




170 


513 




8.89 


687 


0.92 






0.67 


b 


6.1 


380 




185 


576 




12.6 


1116 


0.83 






0.66 


c 


4.1 


446 


381 


123 


598 


559 


7.14 


201 


1.1 


2.4 


0.076 


0.44 


d 


4.5 


471 


412 


132 


620 


585 


8.81 


298 


1.03 


2.3 


0.076 


0.45 



The two first sets (a, b) are obtained in the framework of leading order 
SPA, along lines suggested in our recent work jTJ]. The quantum fluctuations 
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given by the Lagrangian AC r are not taken into account. One can see that 
reasonable fits to the pseudoscalar spectrum at leading order correspond to 
a large value of e ~ 0.7 and clearly show the slow convergence of the 1/N C 
series. 



TABLE 2. The light pseudoscalar nonet characteristics (in units of MeV, except 
for the angle P , which is given in degrees) are presented in a full correspondence 
with the parameter sets of Table 1. 





/- \l/3 

— {uu) ' 


-{ss) 1 


m n 


m K 


u 


fx 


m v 


my 




a 


244 


204 


138* 


494* 


92* 


121* 


487 


958* 


-12.0 


b 


233 


182 


138* 


499* 


92* 


115.8* 


477 


958* 


-15.0 


c 


287* 


263 


138* 


499* 


92* 


115.8 


503* 


958* 


-12.5 


d 


274* 


244 


138* 


494* 


92* 


113* 


494 


958* 


-13.3 



In the sets (c, d) we include the semi-classical correction to the leading 
order result. In these cases a new mass scale parameter A E enters the fitting 
process. Consequently, we have included the value of the light quark con- 
densate as an additional input. One can see that quantum fluctuations lead 
to a small effect, slightly improving the fit. Nevertheless, our expectations 
seem to be realized. The cut-off A E has a large value A E ~ 2.4 GeV, with 
( ~ 0.08 being small enough. 

Let us try to understand why A E — 2.4 GeV and not much larger. The 
reason for this is very simple and is contained in the value for e ~ 0.44. 
This value indicates that cancellations in AC r are not so strong as one would 
expect by using the large- N c arguments to obtain (|19jl. Being 35% less, in 
comparison with sets (a, b), the model parameter e is still too large to vouch 
for suppression in Lagrangian AC r in full measure. 

On the other hand, the calculations would become unreliable if the value 
of A E would lead to such large corrections that they could compete with the 
leading term. We have found that at A E ~ 2.4 GeV the correction amounts 
at most to ~ 20% of the leading term (in the case of weak decay constants 
f n and fx)- This is a strong signal that we are already quite close to exhaust 
the reserves of the quasi-classical expansion and that a further increase of A E 
would start to destroy the fast convergence of the series, giving, in addition, 
a worse fit. 
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4. Discussion of the result 

We have pointed out a new interesting feature related with the 't Hooft deter- 
minant resolution of the Ua{1) problem, which provides one more argument 
in favour of this six-quark interaction: it also unravels the obviously non- 
trivial 7]' puzzle. Indeed, the phenomenological consequences of the 't Hooft 
interaction are well-known. For instance, it leads to spontaneous chiral sym- 
metry breaking with the characteristic mass scale A ~ 1 GeV extracted from 
the quark loops. It also explains the deviations from Zweig's rule. In ac- 
cordance with standard iV c counting rules, the six-quark interaction yields 
a flavour singlet rf meson mass of order m 2 q , ~ 1/N C , as it is commonly ex- 
pected. However, it has not been clear why this 1/N C suppressed mass is 
not much smaller than its actual value of almost one GeV. The mass scale 
for chiral symmetry breaking, A, is too low to explain this phenomenological 
fact. 

To find the answer we have suggested to bosonize the quark determinantal 
interaction and to take into account the next to the lowest order term in the 
semi-classical expansion of the bosonized 't Hooft Lagrangian. This formal 
expansion in powers of the dimensional parameter h actually contains a small 
dimensionless parameter (, which hides the large characteristic scale A E ^ 
2.4 GeV. 

Our solution is obtained in the framework of a simple model, although it 
is strongly based on the instanton picture of the QCD vacuum. Accumulating 
the most essential features of the instanton physics, it is not an accident at 
all, that it resolves the rj' puzzle in a very similar way to the known solution in 
the framework of QCD sum rules. There are, however, important differences: 
in the latter approach the rf gets its mass through mixing with glue states 
and these are expected to be much heavier, giving a new large mass scale 
parameter for the pseudoscalar channel. 

The 't Hooft determinantal interaction is a remnant of gluodynamics at 
large distances, i.e. at scales where quarks interact with each other through 
their zero modes in the instanton background. By means of the semi-classical 
expansion of the bosonized 't Hooft Lagrangian we are able to "touch" the 
border of the non-perturbative region from its low energy side, as opposed to 
the QCD sum rules method. Amusingly, the numerical result is in a perfect 
agreement with our expectations and with the value obtained on the basis of 
QCD sum rules. 



10 



Acknowledgements 



This work has been supported by grants provided by Fundagao para a Ciencia 
e a Tecnologia, POCTI/35304/FIS/2000 and Centro de Ffsica Teorica unit 
535/98. 

References 

[1] S. L. Adler, Phys. Rev. 177 (1969) 2426; J. S. Bell and R. Jackiw, Nuovo 
Cimento 60A (1969) 47; S. L. Adler and W. A. Bardeen, Phys. Rev. 182 
(1969) 1517. 

[2] E. Witten, Nucl. Phys. B156 (1979) 269; G. Veneziano, Nucl. Phys. 
B159 (1979) 213. 

[3] A. M. Polyakov, Phys. Lett. B59 (1975) 82; A. M. Polyakov, Nucl. Phys. 
B120 (1977) 429; A. A. Belavin, A. M. Polyakov, A. Schwartz and Y. 
Tyupkin, Phys. Lett. B59 (1975) 85; G. 't Hooft, Phys. Rev. Lett. 37 
(1976) 8; G. 't Hooft, Phys. Rev. D14 (1976) 3432; Erratum: ibid D18 
(1978) 2199; C. Callan, R. Dashen and D. J. Gross, Phys. Lett. B63 
(1976) 334; R. Jackiw and C. Rebbi, Phys. Rev. Lett. 37 (1976) 172; S. 
Coleman, The uses of instantons, Erice Lectures, 1977. 

[4] D. Diakonov, Chiral symmetry breaking by instantons, Lectures at the 
Enrico Fermi school in Physics, Varenna, June 27 - July 7, 1995, 
-ph/9602375; D. Diakonov, Chiral quark-soliton model, Lectures at 
the Advanced Summer School on Non-Perturbative Field Theory, Penis- 
cola, Spain, June 2-6, 1977, |hep-ph/9802"298| 

[5] G. 't Hooft, Phys. Rev. D14 (1976) 3432; Erratum: ibid D18 (1978) 
2199. 

[6] E. Witten, Nucl. Phys. B160 (1979) 57. 

[7] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; 124 (1961) 
246; V. G. Vaks and A. I. Larkin, Zh. Eksp. Teor. Fiz. 40 (1961) 282. 

[8] P. Di Vecchia and G. Veneziano, Nucl. Phys. B171 (1980) 253; E. Witten, 
Ann. Phys., NY 128 (1980) 363; C. Rosenzweig, J. Schechter and G. 
Trahern, Phys. Rev. D21 (1980) 3388. 



11 



[9] V. Bernard, R. L. Jaffe and U.-G. Meissner, Phys. Lett. B198 (1987) 92. 

[10] V. Bernard, R. L. Jaffe and U.-G. Meissner, Nucl. Phys. B308 (1988) 
753. 

[11] H. Reinhardt and R. Alkofer, Phys. Lett. B207 (1988) 482. 

[12] T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994) 221. 

[13] V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, 
Nucl. Phys. B191 (1981) 301. 



[14] A. A. Osipov and B. Hiller, Phys. Lett. B539 (2002) 76, hep-ph/0204182 



[15] A. A. Os ipov and B. Hiller, Eur. Phys. J. C35 (2004) 223, 
|hep-th/0307035l 

[16] A. A. Osipov, A. H. Blin and B. Hiller, (in preparation). 

[17] A. A. Osipov, H. Hansen and B. Hiller, Nucl. Phys. A (in press), 



hep-ph/0406112 



12 



